Neutralinos, Big Bang Nucleosynthesis and 6 Li in Low-Metallicity Stars 
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The synthesis of 6 Li during the epoch of Big Bang nucleosynthesis (BBN) due to residual anni- 
hilation of dark matter particles is considered. By comparing the predicted 6 Li to observations of 
this isotope in low-metallicity stars, generic constraints on s-wave dark matter annihilation rates 
into quarks, gauge bosons, and Higgs bosons are derived. It may be shown that, for example, 
wino dark matter in anomaly-mediated SUSY breaking scenarios with masses m x < 250 GeV or 
light neutralinos with m x 20 GeV annihilating into light quarks are, taken face value, ruled out. 
These constraints may only be circumvented if significant 6 Li depletion has occurred in all three 
low-metallicity stars in which this isotope has been observed to date. In general, scenarios invoking 
non-thermally generated neutralinos with enhanced annihilation rates for a putative explanation of 
cosmic ray positron or galactic center as well as diffuse background gamma-ray signals by present- 
day neutralino annihilation will have to face a stringent 6 Li overproduction problem. On the other 
hand, it is possible that 6 Li as observed in low-metallicity stars is entirely due to residual dark 
matter annihilation during BBN, even for neutralinos undergoing a standard thermal freeze-out. 



The nature of the ubiquitous cosmological dark mat- 
ter is one of the outstanding questions in cosmology. 
Though there exist a multitude of proposed candidates, 
much consideration has been given to the supersymmet- 
ric (SUSY) neutralino, provided it is the lightest super- 
symmetric particle (LSP). SUSY extensions of the stan- 
dard model are particularly successful in overcoming a 
number of shortfalls of the standard model, such as the 
hierarchy problem and grand unification. Neutralinos are 
appealing since neutral, likely stable, and endowed with 
annihilation/scattering cross sections which (a) makes 
it likely have the right cosmological abundance and (b) 
makes it be detectable in the not-to-distant future by 
either direct- (e.g. scattering in cryogenic detectors) or 
indirect- (e.g. observation of positrons or 7-rays due to 
residual neutralino annihilation in the Galaxy) detection 
means 

In general, the post-freeze-out number-to-radiation en- 
tropy n x /s of a stable and initially abundant particle 
subject to self-annihilation is given by 



Yf = (r 
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where (av) denotes it's thermally averaged annihilation 
rate, H the Hubble expansion rate at freeze-out, n = 0(1) 
for s-wave (p-wave) annihilation, i.e. (av) = cr x~ n with 
x = m x /T, respectively, and it is understood that quan- 
tities are evaluated at the moment of freeze-out (f) it- 
self. Though Eq. could have been obtained by simply 
equating a typical annihilation rate l/Y x dY x /dt with the 
Hubble rate it turns out to be exact. The abundance of 
Eq. Q yields a x contribution to the present critical den- 
sity fl x of 
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where Xf = m x /T at freeze-out, (<7i>)_25 is the annihila- 
tion rate in units of 10~ 25 cm 3 /s at T = m x and g^ H , g$ 



are energy- and entropy- radiation statistical weights at 
freeze-out, respectively. 

Considering freeze-out from annihilations with ther- 
mal equilibrium initial conditions, a particle in the hun- 
dreds of GeV range with (cri>)-25 ~ 0.2 freezing out at 
Xf w 20 when g H « g s w 86.25 due to the known stan- 
dard model degrees of freedom comes close to the re- 
cently by WMAP determined matter density of £l m h 2 ~ 
0.1126+o;°}g} @ (with h the Hubble constant in units of 
100kms _1 Mpc _1 and excluding baryons). In the case of 
mSUGRA models, where the LSP is typically the bino, 
this may be the case in thin strips of parameter space 
when m x < 500 GeV though the bino (av)-25 is typically 
smaller than 0.2 Other parameter space with viable 
neutralino abundances fl m h 2 < 0.113 (corresponding to 
{av} _ 25 > 0.2) in mSUGRA exist in the so-called FOCUS 
point area at several TeV unifying scalar masses mo, with 
the LSP having a large Higgsino component. Leaving the 
fairly restrictive limitation of mSUGRA a wider array of 
possibilities may occur. When gaugino mass unification 
at the GUT scale is dropped [3, the LSP is frequently 
the wino or higgsino, albeit often with an s-wave anni- 
hilation rate of the order (cru)_25 ~ 1 — 100. Similar- 
large annihilation rates may occur in anomaly-mediated 
SUSY breaking schemes (AMSB) 0,110 where the LSP 
is likely the wino, or even in mSUGRA when annihi- 
lation may occur on Higgs poles (e.g. the so-called A- 
funnel region for annihilation on the CP-odd Higgs) for 
m x « TOjy/2 0- It is also possible that the neutralino 
LSP is very light @. 

Many of these cases have in common that for m x not 
too high, thus avoiding fine-tuning for the electroweak 
symmetry breaking to occur, annihilation rates are too 
large to produce £l m h 2 « 0.113 during thermal freeze- 
out. Nevertheless, there is a number of proposed sce- 
narios generating neutralino dark matter either by non- 
thermal means or due to the existence of extra degrees of 
freedom in the early Universe. Q-balls 0, for example, 
non-topological solitons which occur in the spectrum of 



2 



SUSY theories, may be easily formed after an inflation- 
ary epoch of the Universe. As they are made of squarks, 
and are unstable in SUGRA and AMSB, their R-parity 
and baryon-number conserving evaporation at late times 
Td ~ 1 GeV may not only generate many LSP's but also 
the observed cosmic baryon asymmetry. The thus gener- 
ated LSP's may undergo further self-annihilation to reach 
an asymptotic final abundance dependent on Td |10|. 
Similarly, in AMSB scenarios, gravitinos are typically 
heavy ttiq > 50 TeV, such that their decay occurs before 
BBN. In such scenarios, which circumvent the gravitino 
problem, neutralinos are produced in G (and/or mod- 
uli) decays well after the conventional thermal freeze-out 
of neutralinos from equilibrium 0, 0] . When additional 
substantial non-thermal injection of neutralinos at Td 
occurs, their final present day abundance may still be 
calculated via Eq. QJ and Eq. J5J with xt given by 

xj = max[xD,xtf] (3) 

where Xd = tti x /Td, and Xtf quantifies the appropriate 
freeze-out temperature during thermal freeze-out |ll| . 
and where it has been implicitly assumed that the pre- 
annihilation neutralino abundance exceeds that of Eq. (JQ) 
and Eq. (J2J. LSP abundance estimates may also be 
changed (i.e. increased) when compared to the standard 
evolution of the Universe, when additional degrees of free- 
dom are present in the early Universe. This may be, for 
example, due to the existence of a quintessence field in 
the stage of kination, or quite generally, due to extra de- 
grees of freedom gu (with gs unchanged) contributing to 
the Hubble expansion [12| . There are currently no limits 
on gn at high T other than those from the (much later) 
epoch at BBN T < 1 MeV. 

Though the post-freeze-out Y£ stays essentially un- 
changed, residual x annihilations occur up to the present 
day and may have considerable impact on the process 
of Big Bang nucleosynthesis (BBN). The effects of resid- 
ual annihilation ofy's during BBN have been considered 
some time ago [lj, [l4|. Nevertheless, whereas earlier 
studies were based on the isotopes 2 H and 4 He , it has 
been poin ted out that 6 Li production during electromag- 
netic [liJ Ha] or hadronic [l^, [13, Ha energy injection 
during BBN is particularly efficient [2(j. The isotope of 
6 Li may be produced via non-thermal reactions (subject 
to energy threshold), such as 3 H(a, n) 6 Li or 3 He(a,p) 6 Li. 
Energetic 3 H and 3 He are readily produced via nuclear 
4 He spallation by energetic nucleons or 4 He photodisin- 
tegration. Since the photodisintegration of 4 He is only 
efficient for T< 0.3 keV particular importance in the an- 
nihilation case is due to hadronic spallation. This is be- 
cause residual annihilation is stronger at earlier times. 
Resulting 6 Li abundances are given formally by a time 
integral 
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where the factors on the right-hand-side from left to 



right, denote (a) the x annihilation rate 

dAnn _ (<™> f Px V fK\ 

with p x neutralino density, (b) the generated energetic 
protons p and neutrons n per annihilation, (c) the pro- 
duced energetic mass-three nuclei per generated ener- 
getic nucleon, and (d) the produced 6 Li per energetic 
mass-three nucleus, respectively. Here the first factor is 
determined by the annihilation rate under the assump- 
tion that x is the dark matter whereas the second fac- 
tor may be obtained assuming a particular x annihila- 
tion channel and computing the energetic nucleon spec- 
trum via a hadronic flux tube Monte Carlo code such 
as PYTHIA [2l|. Concerning the third and fourth fac- 
tors I have recently presented first ^3 results and a first 
description of a newly developed Monte-Carlo code de- 
scribing the cascading of energetic nucleons and nuclei 
on background thermal nucleons, nuclei, and electromag- 
netically interacting particles which allows me to evalu- 
ate Eq. (0J. This code includes, Coulomb and Thom- 
son stopping of fast charged nucleons and nuclei, elastic- 
and inelastic- nucleon- nucleon scattering and elastic- and 
break-up- processes in nucleon- 4 He scattering. Of all 
mass three nuclei produced a small fraction react on 4 He 
to form 6 Li. The probability to do so, may be approxi- 
mated by 

P ^ = (-) = [ dE ^iw f dE i 4 

\n 3 J J &E t J Eth U u {dE/dx) c 

xexpf- f *dE'- * ) (6) 

V Je lnuc(dE/dx)c J 

a convolution over the initial energy distribution of mass- 
three nuclei, dn^/dEi, - the probability that a mass-three 
nucleus during it's passage through energy space due to 
the dominant multiple Coulomb interactions with energy 
loss per unit length dE/dx undergoes a reaction to form 
Li (where Ze Li = l/(c76Lin4 He ) is the mean free path to- 
wards formation of 6 Li) - while, the exponential factor, 
not having undergone already another nuclear interac- 
tion (mean free path due to all nuclear interactions is 
Inuc) 10] ■ Compared to the Coulomb losses other nu- 
clear interactions in Eq. © are typically not very impor- 
tant, except for charge exchange between the mirror nu- 
clei 3 H(p, n) 3 He and to a lesser degree elastic p- 3 H and 
p- 3 He scattering, in only a narrow temperature range 
between several keV and « 20 keV. Nevertheless, though 
the freshly synthesized energetic 6 Li typically survives 
almost completely p-spallation p3| it may only survive 
6 Li(p, a) 3 He thermal reactions (after it's rapid thermal- 
ization) when T< lOkeV. It has been realized [T^L ITflj 
that in the above-given temperature window Coulomb 
losses loose some of their efficiency (by roughly a factor 
3-10) for energetic 3-nuclei with velocities below the elec- 
tron thermal velocities (but still above the lithium forma- 
tion reaction thresholds of 8.39 and 7.05 MeV for 3 H and 
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FIG. 1: Final 6 Li yield functions defined via Eq. as a 
function of neutralino mass for various annihilation channels 
as labeled in the key. The 1-a range of the 6 Li abundance in 
HD84937 is also shown. 

He, respectively). I have computed the energy transfer 
of a fast charged nucleon or nucleus due to Coulomb in- 
teractions with an electron-positron plasma at tempera- 
ture T, explicitly accounting for a thermal average and 
accounting for higher order terms. Details on this cal- 
culation will be presented elsewhere. Since 6 Li synthe- 
sis during/after BBN is dominated during those epoch 
where (a) the annihilation rate Eq. JSJ is still large, (b) 
synthesized 6 Li survives the thermal 6 Li(p, a) 3 He reac- 
tion, and (c) dE/dx\c is at it's minimum, the bulk of 
the 6 Li is synthesized at T « 10 keV and an exact eval- 
uation of dE/dx\c is paramount to an evaluation of 6 Li 
abundances. 

I have thus computed the 6 Li yield for annihilating par- 
ticles and under the assumption of specific annihilation 
channels such as into uu, dd, ss, bb, or ti quark-antiquark 
pairs, as well as W~W + or ZZ gauge bosons [24| ■ The re- 
sults of these computations for varying neutralino masses 
m x are shown in Fig. 1. Here the annihilation-channel 
dependent yield functions Y e l Li arc defined via the equa- 
tion 

n6 U , \ -3/2 ( Q X h2 \ S^uvi (>7\ 

— = M-25 m x , 100 ^_ j J2 W u (7) 

where the bi are branching ratios into channel i and 
m x,wo is the neutralino mass in units of 100 GeV. Eq. Q 
is remarkable as it allows a quite general evaluation of 
the final 6 Li abundance, independent of the nature of 
the annihilating particle, as well as applicable for essen- 
tially all relevant (av) and m x . Note that the Ye Li 's for 
leptonic channels, which are not shown, are essentially 
zero due to the absence of injected nucleons. Annihi- 
lation into Higgs bosons, on the other hand, yield Ye Li 
only somewhat smaller than those shown in Fig.l, since 
Higgs bosons typically decay into heavy quarks or mas- 



io- 23 


3 

<0"v> (cm /s) 


T D = 100 MeV ' 




8g/g = 10 4 




io- 24 


5g/g= 10 2 




10 




T D = 5 0eV 


IO" 26 


y"^ standard thermal freeze-out 

- 

M x (GeV) 




10 


100 1000 



FIG. 2: Annihilation-channel dependent constraints on the 
s-wave annihilation rate due to possible 6 Li overproduction 
as a function of neutralino mass. For simplicity, only the uu, 
bb and W~W + channels are shown with heavy lines and line 
styles as indicated in Fig.l. Constraints on other quark- and 
gauge boson- annihilation channels are virtually identical to 
those shown as may be verified by inspection of Fig.l. Anni- 
hilation rates above the lines are ruled out. Also shown are 
the annihilation rates required [3] to produce Q, x h 2 = 0.1126 
during standard thermal freeze-out (solid line) or during ther- 
mal freeze-out when extra degrees of freedom (contributing 
to the Hubble expansion during freeze-out) are present (the 
upper two solid lines with Sg^/g^ as labeled). The dotted di- 
agonal lines correspond to the required (av) for post thermal 
freeze-out non-thermal generation of Q, x h 2 = 0.1126 at tem- 
peratures, from top to bottom, 0.1,0.2,0.5,1,2, and 5 GeV, 
respectively. Here the QCD phase transition has been as- 
sumed to occur at 200 MeV. The curved dotted line shows 
the annihilation rate in case of AMSB winos |5j. 



sive gauge bosons. A simple scaling with (av) was possi- 
ble due to the linearity in the 6 Li production and destruc- 
tion mechanisms. Note that though for fixed dark matter 
density d Ann/di ~ m~ 2 a scaling in Eq. J7J) with m^ 15 
has been adopted. This is due to higher mass x P r0_ 
ducing more energetic primary nucleons which, in turn, 
produce a larger number of secondary p and n, yielding 
a final ne Li /?i,H ~ m~ 15 . 

Fig. 1 also shows the one-sigma range 6 Li/H« 8.47 ± 
3.10 x 10~ 12 j2|| of the observationally best determined 
6 Li/H-ratio in a low-metallicity star, i.e. in HD84937, 
which has been analyzed by several groups jlEESl- 6 Li 
detections have been currently claimed in three low- 
metallicity [Z] < — 2 12a. |28l |29| and two higher metal- 
licity stars [Z] ~ —0.6 30] with their abundances coin- 
cidentally all in the same range, reminiscent of a unique 
cosmic abundance. In principle, 6 Li may be depleted 
in stars, in practice, however, those low [Z] stars which 
show 6 Li have normal "Spite" -plateau 7 Li abundances, 
and as this latter isotope would, in most circumstances, 
be depleted as well, substantial 6 Li depletion seems un- 
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FIG. 3: S-wave annihilation rate required to produce within 
the 2 — o limits the 6 Li abundance of HD84937. The heavy 
lines indicate the central value of HD84937, whereas lighter 
lines the 2 — a ranges. For simplicity only the uu (solid) and 
W~W + (dotted) channels are shown with results for other 
channels similar (cf. Fig. 1). 

likely E]. The origin of the 6 Li at [Z] < — 2 is somewhat 
mysterious as only with considerable difficulty explained 
by traditional galactic cosmic ray spallation and fusion 
reactions It is thus conceivable that 6 Li at low [Z] 
is, in fact, entirely of primordial origin, though other al- 
ternative origins have also been proposed [33j . 

In the absence of stellar 6 Li depletion the efficient pro- 
duction of 6 Li due to annihilating neutralinos may be 
used to constrain the properties of the dark matter par- 
ticle itself. Fig. 2 shows annihilation-channel dependent 
limits on the s-wave annihilation \X cross section of a 
dark matter particle, representing a convolution of the 
results in Fig. I with Eq. (JJJ. Here the two-sigma up- 
per limit on 6 Li in HD84937 has been adopted such that 
dark matter particles above the diagonal lines in the up- 
per part of the figure are ruled out due to 6 Li overpro- 
duction. For reference, the figure also shows annihilation 
cross sections required to obtain the WMAP abundance 
given (a) standard thermal freeze-out 01 (b) thermal 
freeze-out when additional degrees of freedom are present 
5g H /g H > as, for example due to a quintessence field 
and (c) via late generation below the thermal freeze-out 
temperature due to, for example, decay of gravitinos or 
evaporation of Q-balls. Furthermore, the figure shows 
the predicted LSP wino >(av) into W + W~ in AMSB sce- 
narios as given in Ref. |5j 

It is evident that possible 6 Li overproduction imposes 
stringent constraints on dark matter particle proper- 
ties. As evident from Fig. 2, light neutralinos with mass 
m x < 20 GeV annihilating into light quarks are ruled out 
at two-sigma. AMSB winos may only be consistent with 
the 6 Li abundance when m x > 250 GeV. When the dark 
matter particle is the first excited Kaluza-Klein mode 



B^ j3|| of the U(1)y gauge boson, it may not be 
lighter than m^K $ 200 GeV due to its efficient annihi- 
lation into right-handed up- type quarks. Late-time gen- 
eration of dark matter x's is disallowed for all tempera- 
tures below 200 MeV unless the dark matter particle is 
fairly heavy m x > 400 GeV or annihilation does not occur 
into hadronic or gauge (and higgs) boson channels. Fi- 
nally, a substantially increased Hubble rate during ther- 
mal freeze-out is for all but the heaviest neutralinos ob- 
servationally disallowed. These limits may be circum- 
vented in case substantial stellar 6 Li depletion occurred 
which, nevertheless, would require a coincidentally sim- 
ilar amount of depletion in all three 6 Li-rich observed 
stars to date. Limits of this sort are also important in 
light of scenarios which invoke neutralino annihilation 
as putative explanations of, for example, the observed 
cosmic positron excess at ~ 10 — 30 GeV I37ll38t a s de- 
termined by HEAT, or the galactic-center \m 14(1 l4l| . 
or extragalactic diffuse J42j , gamma-ray fluxes as deter- 
mined by EGRET, VERITAS, or CANGAROO. In order 
to explain anomalous components of such signals such 
as bumps in the spectrum a signal boost (enhancement) 
factor B s of the order - 50 - 10 3 [H Ej is essentially 
always required. Considering the most recent N-body 
simulations on substructures and halo profiles ^fjj such 
B s seems only unlikely due to clumpy halos or singular 
halo profiles though, in principle, one could envision it 
due to an enhanced (av) with respect to it's standard 
thermal freeze-out value. Nevertheless, already modest 
particle-physics motivated boost-factors of the order of 
~ 1 — 10 will have to face a potential 6 Li overproduction 
problem. 

Last but not least, it is possible that the entire ob- 
served 6 Li at low metallicity may be due to the residual 
annihilation of a dark matter particle. Fig. 3 shows the 
mass-dependent annihilation rate required to produce a 
6 Li abundance within the 2a ranges of those observed 
in HD84937. It is seen that this may be accomplished 
even by a standard thermal freeze-out with dominant s- 
wave component annihilation into light quarks, provided 
the neutralino mass is within the approximate range of 
20 - 80 GeV 0. The observed amount of 6 Li may be 
produced for even larger mass neutralinos when either 
coannihilation effects or annihilation on poles occur in 
the thermal freeze-out case or neutralinos are generated 
non-thermally. Coincidentally, the recently proposed 
specific dark matter neutralinos 0] and Kaluza-Klein 
particles |4?| which could explain the claimed bump in 
the extra-galactic 7-ray background and/or the positron 
excess as observed by HEAT would have just the right 
properties to yield 6 Li abundances as observed in low- 
[Z] stars. Such annihilation is also associated with some, 
albeit small, amount of observationally favored 7 Li deple- 
tion 

El- 

It is intriguing that the observed abundances 
of 6 Li in low-metallicity stars may be entirely a product 
of dark matter annihilation. 

I acknowledge discussions with Eric Nuss. 
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